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^ ■ Abstract 

, The equation for the Wigner function describing the reduced dynamics of a single harmonic 

O 

^ \ oscillator, coupled to an oscillator bath, was obtained by Karrlein and Grabert [Phys. Rev. E 55, 
153 (1997)]. It was shown that for some special correlated initial conditions the equation reduces, 

■ in the classical limit, to the corresponding classical Fokker-Planck equation obtained by Adelman 

, [J. Chem Phys. 64, 124 (1976)]. However for separable initial conditions the Adelman equations 

g 

I ' were not recovered. We resolve this problem by showing that, for separable initial conditions, the 

^ ■ classical Langevin equation obtained from the oscillator bath model is somewhat different from the 

o 

o : one considered by Adelman. We obtain the corresponding Fokker-Planck equation and show that 

it exactly matches the classical limit of the equation for the Wigner function obtained from the 

> 

lO ' master equation for separable initial conditions. We also discuss why the special correlated initial 

^ : 

. conditions correspond to Adelman's solution. 

I ■ 
I ■ 
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I. INTRODUCTION 



The concept of 'open' quantum systems is a ubiquitous one in that any real system is 
never fully isolated. It is usually in contact with a larger system, constituting its environ- 
ment, which effects its dynamics. Ford, Kac and MazurA suggested the first microscopic 
model describing dissipative effects. They considered a system which was coupled to a bath 
modelled by a collection of an infinite number of harmonic oscillators. Interest in quantum 
dissipation was intensified by the work of Caldeira and Legget^ who used the influence func- 
tional approach, developed by Feynman and Vernon^, to discuss quantum Brownian motion. 
They considered the case where the system and its environment were initially uncorrelated, 
the so called separable initial condition. This was generalized to the situation where initial 
correlations exist between the system and its environment by Hakim and Ambegaokai^, 
Smith and Caldeira^, Grabert, Schramm and Ingold®, Chen, Lebowitz and Liveranii and 
Banerjee and Ghosh^ among others. The master equation for the quantum Brownian motion 
of a harmonic oscillator in a bosonic bath of harmonic oscillators was obtained by Haake 
and Reibold^ and later by Hu, Paz and Zhangi^ who used path integral methods in their 
derivation. A derivation of the Wigner transform of this equation was also given by Halliwell 
and Yi>ii and a solution of this Fokker-Planck equation with time dependent coefficients was 
given by Ford and O'Connell-'^-. 

In a very interesting paper Karrlein and Graberlii^ studied the master equation for the 
reduced density matrix of the system for generalized correlated initial conditions. They 
showed that in general it is not possible to obtain a Liouvillean operator for the master 
equation of the reduced dynamics of the system under the influence of its environment. 
However for some specialized correlated initial conditions like the one discussed by Hakim 
and Ambegaokari, called the thermal initial condition, and for the case of the time evolution 
of equilibrium correlation functions®, it is possible to obtain the Liouvillean to describe the 
reduced system dynamics. In both these cases Karrlein and Graberlji^ obtained the Fokker- 
Planck equation from the Wigner transform of the master equation. They found that in 
the classical limit, these equations reduced exactly to the corresponding equation obtained 
by Adelmani^ for the dissipative dynamics of a classical harmonic oscillator driven by a 
Gaussian non-Markovian noise. They then considered the case of separable initial conditions 
and found that in the classical limit the corresponding Fokker-Planck equation did not reduce 



2 



to the Adelman equation. We find this a very intriguing as well as an interesting point and 
the main purpose of this paper is to look at this point more closely. This is also important 
to understand in view of the fact that factorizing initial conditions are known to have some 
problematic features-^^ associated with the sudden switching on of the coupling. 

The plan of this paper is as follows. In Sec. II we obtain the Langevin equation for the 
system consisting of a harmonic oscillator coupled to a bath of harmonic oscillators with 
which it is initially uncorrelated. We indicate that our treatment is valid in the classical 
as well as the quantum regime. In Sec. Ill we obtain the Fokker-Planck equation from the 
Langevin equation derived in Sec. II in the classical regime. When we compare our equation 
with the one used by Adelman-- we find that our equation has an additional term. This 
additional term is the crux which enables us to resolve the issue. The Fokker-Planck equation 
obtained from the Langevin equation with this additional term is the exact classical limit 
of the corresponding equation in the quantum regime obtained by Karrlein and Graberlii^. 
In Sec. IV we discuss why the Langevin equation used by Adelman would be the correct 
classical equation describing a correlated initial state. In Sec. V we discuss our conclusions. 



II. LANGEVIN EQUATION 

We consider a single harmonic oscillator described by the phase space variables {x,p) 
and with a natural frequncy ujq. The oscillator is coupled to a heat bath modelled by 

independent oscillators described by {X^, Pa} and which have frequencies {00^} ■ To 
generate dissipation one eventually takes the limit A^ ^ 00. The full Hamiltonian is taken 
to be: 

H = '^ + ^ + U^ + ^]-j:c.xXa. (1) 

^ ^ a=l ^ ^ a=l 

The last term represents the coupling between the system and the bath. Introducing the 
shifted frequency u"^ = Uq — Z^qC^/u;^, we can rewrite Eq. (H)) in the following equivalent 
form 



2 2 2 N p2 ujliXa-^f 
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Let us now write the equations of motion for the system and bath variables. These are given 
by: 

X = -ulx + Cc^Xc, , (3) 

a 

Xa = -UjlXa + CaX . (4) 

The solution of the equations of motion of the reservoir variables gives 

Xa{t) = cos(cj„t)X„(0) + — ^^X„(0) + / dt' ^ -Cax{t') . (5) 

Ua JO ijJa 

Plugging these into the system's equation of motion Eq. Q we get the following Langevin 
equations of motion: 



X = -ulx + [ dt'^^^^y^xit') + r^it) , (6) 



^2 







where 7(t) = XI ~| cos^Uat) , (7) 

Vit) = Eca[cosKt)X„(0) + ^^^i^^^X„(0)] . (8) 

Note that the discussion till now is valid both clasically and quantum mechanically. In 
the quantum case the variables correspond to operators in the Heisenberg representation. 
We now use the fact that at t = the separable initial conditions imply that the reservoir 
(isolated) is initially in equilibrium and hence we know the statistical properties of the bath 
variables {Xq,(0), Xa(0)}. In the quantum case we get 

2v2/'n\ \ I <r2ir.\ \ 'huJa 



(<^.(o)) = (x^(o)) = Y^oth.^^^^, 

( [ X,(0)X„(0) + X„(0)X,(0) ] ) = 0. (9) 

In the classical case we get, using equipartition, (tu^X^) = (X^) = ksT and (X^Xq) = 
and these also follow from Eq. Q if we take the classical limit huOa/ (ksT) — * 0. From these 
properties of the bath one can work out the correlations of the noise ri{t) in Eq. (jH}. We get 
(for t, t' > 0) 

( vit) ) = , 

li vitMt') + vit'Mt) ) = hK'it-t'), (10) 
where K'it) = ^ A cos coth . (11) 



Here we have used the notation in Ref . |6ljl3l| . The classical limit of Eq. ()1U|) gives the usual 
fluctuation-dissipation relation 



(12) 



III. FOKKER-PLANCK EQUATION FOR THE CLASSICAL LANGEVIN EQUA- 
TION 

In this section we focus on the classical Langevin equation which is described by Eqs. 
First we note that, after a partial integration, we can rewrite Eq. (jU)) in the following 

form: 

X = -u^x -^{t)x{0) - dt'^{t -t')x{t') +r]{t) . (13) 

This differs from the equation considered by Adelman^'* by the initial condition dependent 
term — 7(t)a;(0), the so called initial slip ter m^^i^^ which, as we will see, is crucial in get- 
ting the correct Fokker-Planck equation. Starting from the above Langevin equation or 
its equivalent form Eq. © we get, following Adelman, the corresponding Fokker-Planck 
equation. 

The general solution of Eq. (jU)) is given by: 

x{t) = H{t)x{0) + G{t)p{0) + f dt'G{t - t')T]{t') , (14) 

Jo 

where H{t) and G{t) are solutions of the homogeneous part of Eq. (jHI) for the initial condi- 
tions H{0) = 1, if (0) = and G(0) = 0, G{0) = 1 respectively, and p{t) = x{t). Taking a 
Laplace transform of the homogeneous equation we immediately get: 

^ 1~, — ' 

^(^) = ~, — ^ ' 

Z"^ + 00"^ + Z ^[Z) 

where G{z), H{z) and 7(z) are the Laplace transforms of G{t), H{t) and 7(t) respectively. 
We note that G{t) is the same Green's function, denoted by Xuif), in Adelman'si^ paper. 
However H{t) is different from Adelman's Xx) which we will here denote by -F(t), and which 
is given by 

F{t) = 1-00^ f' dt'G{t') . (17) 

JO 
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From the form of the Laplace transforms G, H and from Eq. (jl7|) it is also easy to get the 
relations 

G{t) = Hit) and F{t) = - u^G{t) . (18) 

From the general solution for the position we get for the momentum 

p{t) = H{t)x{0) + G{t)p{0) + r dt'G{t - t')7]{t') . (19) 

Jo 

Let us use the notation {xi, X2} = {x, p}. The fluctuations about the mean position and 
momentum are given by: 

yi{t) = xi{t) - [H{t)x{0) + G{t)pm = f dt'G{t-t')r]{t') , 

Jo 

y^it) = X2{t) - [H{t)x{Q) + G'(t)p(0)] = f dt'G{t - t')r]{t') . (20) 

Jo 

In the classical case, since ri{t) depends linearly on Xa{0) and Xa{0) which have a Gaussian 
distribution, it follows that ?7(t) is also Gaussian. Hence the distribution of {^1,1/2} at any 
time will also be Gaussian and is completely specified by the following variances, obtained 
from Eqs. (fT ^ . 



Au{t) = {yf) = ksT t dt' f dt" G{t') G{t") 7(t' - t") 

Jo Jo 

A22{t) = {yl) = keT fdt' fdt" Git') Git") 7(t' - t") 



Jo 



Ai2it) = A2iit) = {yi y2) = ksT f dt' f dt" Git') Git") 7(t' - 1") . (21) 

JO Jo 

Taking time-derivatives and after some manipulations one can get the following results for 
the variances: 

FF . 



An = -2kBT iGG + 

F^-1 



Cj2 



An = -ksT i G 



2 



Cj2 



A22 = -2kBT iG G + uj^ GG) , 
A22 = -ksT i uj^G^ + G^-1) , 

= keT iG F-GG) . (22) 

These are of course identical to the results of Adelman since the variances involve exactly 
the same Green's function. However the mean values (a;(t)), {pit)) are different and we will 



now show that this leads to a different Fokker- Planck equation. The phase space probabihty 
distribution function for {x, p} is thus 

P{x,p,t\x{0),piO),0) = 2n{det[A]y/- '~"^^'"'^' " ^^^^ 
Following Adelman, let us define the following matrix: 

/ 



H{t) G{t) 
m G{t) 



(24) 



and new variables {qi, ^2} through the transformation 



q2it)) \X2{t) 

Clearly we then have qi{0) = x(0), ^2(0) = p(0). In terms of these new variables the 
probability density P takes the following form: 

P(gi,g2,tk(0),p(0),0) = P{x,,X2,t\x{0),p{0),0) 

1 



27r(det[A])V2 
where B = TAT^ . 



te-«i(o)]B-.'te-%(o)]/2 ^ ^26) 



Taking time and space derivatives of P gives the following relations: 
BP 1 (7 _ 1 . 

l^^^irn = —BfJiBi^P + ^Bf,' B-] Bi^ [q, - g,(0)] - g,(0)] P 

= I In det[B] P - 1 Bf^ [q^ - ^,(0)] [g^ - g^(0)] P , (28) 
where the last line follows from using the following identities for the symmetric matrix B: 

-det[B] = ^ ^ Bim = det[B] Bf^ Bi^ , 

at 0-Dlm 

B"^ = -B^BB"^ . 

We will also need the following result: 

dP _ dP • _i dP 
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Using Eqs. (j2ZI), (PHj) . (1221) transforming from the variables {qi ,^2} to {xi, X2} we 
finally get: 



f ^-|t-Tl,„..„|;-|.„AP + l|T-BT-- 



\lm 



dxidxr 



(30) 



where A = det[T-^] 
and get: 

T 1 t = 



— G G. We work out the explicit forms of the various matrices 



1 







where 7^ 

7p 



A + T"^ t A + (T-^ t A)^ 

^ in - 2Ai2 ii2 - ^22 + Iq All + 7p M2 

^ il2 - ^22 + Iq All + 7p ^12 ^22 + 27gAi2 + 27p ^22 



G2 - G G 
GG-G G 



— In A . 

dt 



G^-G G 

From Eq. we get An = 2 Ayi- Using this and expanding out Eq. (jHUjl gives: 



(31) 
(32) 

(33) 

(34) 
(35) 



aP(a;,p,t|x(0),p(0),0) 



dP 



dP d{pP) 



-PjT + 7g a; ^ + 7p- ^ 
ox op op 



d^p 

dp^ 



d^p 

dxdp 



(36) 



1 ■■ 7p • 

where Dq = -An - A22 + 7g ^11 + y ^11 



1 • 7 • 

-A22 + -J An + Ip A22 



From Eq. (|7j) and Eq. (fTTj) it is clear that in the classical limit: 

K\t) - ^7(t) , 



(37) 



which implies that the classical limits of the functions Kq{t) and Kp{t) in Ref. 13| is: 



h Kq{t) An, h Kp{t) ^ A22 . 



(38) 



With this we immediately see that the coefficients of the equation for the Wigner function, 
given by Eq. (89) in Ref. Q], reduces in the classical limit, to those in the Fokker-Planck 
equation given by Eq. (jHUj) . 
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IV. CORRELATED EQUILIBRIUM INITIAL CONDITIONS 



It is possible to understand the precise reason for the observed correspondence between 
the Adelman results and the quantum Fokker-Planck equation for the correlated initial 
conditions considered in ^|. We will now show that the Langevin equation used by Adelman 
also follows from the same microscopic model if, instead of choosing a separable initial 
condition, we take a correlated initial state which is prepared by making measurements on 
an equilibrium state of the coupled system and bath. We first make some comments on the 
physical significance of the special correlated initial conditions that have been studied in 
the quantum case. The two initial conditions correspond to (i) making measurements of the 
position variable of the system at time t = on the equilibrium state, (ii) computing time 
dependent correlation functions of various system variables, again in an equilibrium state. 
Especially it is to be noted that the measurements are made on the equilibrium state of the 
fully coupled system and hath. Classically if we want to describe a similar situation, then the 
corresponding Langevin equation should decribe the time evolution of the system's phase 
space variables given that an initial measurement, in the initial equilibrium state of 

the coupled system and bath, gave the values {x(0),p(0)}. 

Starting from the same model of reservoir with the full Hamiltonian given by Eq. or 
equivalently Eq. (0) we again get an equation in the form of Eq. (fT^ but the noise correlations 
are not given by Eq. ()12j). At t = the bath and system variables are correlated, hence 
( x(0) Xq,(0) ) 7^ and the bath correlations are no longer give by Eq. (jH]) or its classical 
counterpart. A simplification now occurs if we redefine the noise^^ and write Eq. (|13|) in the 
following form: 

X = -u'^x- f dt'-f{t-t')x{t') + ^{t) , (39) 
where i{t) = r]{t) - 7(t)x(0) 

= Yl Ca{ cosKt) [ X„(0) - ^ x(0) ] + ^l^^i^^X,(0) } . (40) 

The statistical properties of the noise should be evaluated by averaging over the bath vari- 
ables. Following the discussion in the previous paragraph it is clear that the distribution of 
the bath variables is given by the conditional probability of {Xq,(0), Pq,(0)} given that an ini- 
tial measurement, in the initial equilibrium state, on the system gave the values {x{0),p{0)}. 
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From the form of the Hamiltonian Eq. we get for the conditional distribution: 

Pr({X,(0),P,(0)}) 



Zb 



where Hb = Y.I + 



,2 2 

a 

Averaging the bath variables over this distribution we then get: 

([X,(0)-^x(0)]X,(0)) = 0, 

( ^l[Xa{0) - ^ X{0) r ) = kBT, 

( Xl{Q) ) = kBT. (41) 

Using these we then find that ^{t) is again a Gaussian stationary process with exactly the 
same correlations as obtained for ri{t) for separable initial conditions. Thus 

m) = , 

mm) = kBT^it-t') . (42) 

The Langevin equation given by Eq. (jH^ with the noise given by Eq. ()42|) is the starting 
equation of Adelman^. Hence we have shown that the Langevin equation of Adelman is 
the correct classical equation desribing a correlated initial state. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have taken up the issue of the classical limit of the quantum Fokker- 
Planck equation for a harmonic oscillator linearly coupled to a bath of harmonic oscillators 
wherein the system and bath are initially uncorrelated. Karrlein and Graberlii^ obtained 
the quantum Fokker-Planck equation for some special correlated as well as separable initial 
conditions. Two kinds of correlated initial conditions were studied, the so called thermal 
initial condition*^ and those describing the dynamics of equilibrium correlation functions^. 
On taking the classical limit of these equations they found that they exactly reduced to 
Adelman's classical Fokker-Planck equatioriF^. However for separable initial conditions they 
found that the classical limit of their equation did not give Adelman's equation. It was 
suggested that this indicated a problem with separable initial conditions. 
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We have found a simple resolution to this somewhat intriguing point. We derived the 
classical Langevin equations of motion corresponding to the system of a single oscillator 
coupled to a bath of harmonic oscillators and with separable initial conditions. This equation 
differs from the one used by Adelman-^* by an extra piece —'~f{t)x{0) which depends on initial 
conditions. We then derived the Fokker-Planck equation for this Langevin equation and 
found it to be the exact classical limit of the quantum equation of Karrlein and GraberU^ 
for separable initial conditions. We have also shown why for the case of correlated initial 
conditions the Wigner function equation of the reduced system dynamics in the quantum 
regimeiSi reduces to the Adelman equation in the classical limit. We have shown that the 
Langevin equation used by Adelman corresponds to the situation where the initial state 
is the equilibrium state of the coupled system and reservoir. Thus the basic point is that 
correlated and uncorrelated initial conditions lead to different classical Langevin equations 
and hence to different Fokker-Planck equations. In conclusion we have shown that for 
the case of a single harmonic oscillator coupled to a bath of harmonic oscillators there is 
an exact correspondence between the classical Fokker-Planck equation and the quantum 
equations for the Wigner function both for some special correlated as well as uncorrelated 
initial conditions. 
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